Phase synchronization from noisy univariate signals 
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We present methods for detecting phase synchronization of two unidirectionally coupled, self- 
sustained noisy oscillators from a signal of the driven oscillator alone. One method detects soft, 
another hard phase locking. Both are applied to the problem of detecting phase synchronization in 
von Karman vortex flow meters. 
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Phase synchronization of nonlinear, self-sustained os- 
cillators is known since the observations of Huygens over 
300 years ago 0. For several reasons, interest in the phe- 
nomenon has recently resurged, for review, see |^ ^, 0. 
For chaotic oscillators, phase synchronization was found 
to be an independent regime among various others with 
different degrees of synchronization jj^, || [?], Phase 
synchronization has also been observed in several natu- 
ral phenomena [0], most notably perhaps in brain physi- 
ology H [l0[ , where it is thought to play a key role in 
information processing Therefore, interest has risen 
in the problem of identifying and characterizing phase 
synchronization from measured time series [ fl2| . 

For the situation that signals from both oscillators 
are available, several techniques have been developed to 
quantify the degree of phase synchronization between 
them H [n|, to determine the coupling direction 14, ]1| 
and the type of coupling Q. 

For univariate signals from bivariate systems it is 
sometimes possible to identify and separate oscillatory 
components of the two oscillators, and to proceed as in 
the bivariate case JT|, [l?], [ll| . 

But when a univariate signal contains oscillations from 
a single oscillator only, the question whether it is phase- 
synchronized to some other oscillator is much harder to 
answer. Important progress towards this goal is the 
method of Janson et. al. [[l9), where the driving oscil- 
lator is identified by observing the regular variations it 
induces in the period of the driven oscillator. Effectively, 
these observations are used to reconstruct the dynamics 
of the pair of coupled oscillators in the spirit of Takens' 
delay embedding. As for all embedding techniques, the 
method is sensitive to noise, and careful pre-processing 
of the data is required for noise reduction. 

In this Letter, a different method for detecting syn- 
chronization in terms of phase locking from univariate 
time series, based on nonequilibrium phase statistics, is 
proposed. It applies to situations with unidirectional 
coupling where only a signal from the driven oscillator 
is available. Similar to the method of Rosenblum and 
Pikovsky for the identification of the coupling direction 
juj, it is not only robust to noise, but requires the dy- 
namics of the driven oscillator to be perturbed by weak 
internal noise. It turns out that a regime of "hard" phase 



locking has to be distinguished from a regime of "soft" 
phase locking, and two different tests have to be used. 

An advantage of the tests proposed here is that they 
cover all frequency ratios u>\ : loi ~ n : m of driven 
Oscillator 1 to driving Oscillator 2 with small integers n 
and m. This includes the cases with n = 1, where the 
method described in Jl9| cannot be applied. In particular 
for the case of 1 : 1 phase locking, where a separation of 
the signals of the oscillators involved is excluded, no test 
seems to be known so far. 

For simplicity, the theory is first described for this im- 
portant case only. We specify the conditions for a suc- 
cessful detection of phase locking, and apply it to simu- 
lated data. Finally, a demonstration of the method for 
vortex flow meters is reported. 

In most experiments not the phase itself but some oscil- 
latory signal x(t) is measured. One possibility to extract 
the phase <f>(t) from such a signal is to first calculate the 
Hilbert transform xh of the signal, 



XH {t)= i - P .v. r ^idt 



where P.V. denotes the Cauchy principal value, 
the analytical signal 

C(i) = x(t) + ix H (t) = A(t) exp(z0(i)) 



(1) 



Then, 



(2) 



gives a consistent way to define the phase <f>(t) |3[ [l8| . By 
extracting the phase <f>{t), the higher dimensional dynam- 
ics of an oscillating system is projected to one dimension. 

A minimal model || that generally describes the phase 
dynamics of two unidirectionally coupled, self-sustained 
oscillators is given by: 
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-Ksin(A0) +£i 
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where £i and £2 represent Gaussian noise with 
(*0) = 2Dk5k,i5(t - t'), k, I = 1,2 



(3a) 
(3b) 

(3c) 



Unidirectional coupling of two self-sustained oscillators in 
the presence of noise is controlled essentially by the four 
characteristic time scales of the phase dynamics, which 
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FIG. 1: Typical, mean drift corrected, phase 4>i{i) — (d(j>i/dt) t 
simulated with equation || for (a) no, (b) hard and (c) soft 
phase locking and the distributions of the phase increment 
ip T0 for each case (d)-(f). (g) the distribution u{(3), cf. Eq. (0), 
corresponding to (f). 

correspond to the four rate constants entering the min- 
imal model. These are the detuning Aw = uji — u> 2 be- 
tween the oscillators, the diffusion coefficients D\ and 
D 2 of the phases (f>\ and <j> 2 of the two oscillators in the 
uncoupled limit, and the relaxation rate K of the rela- 
tive phase A<p(t) := 4>i(t) — 02 (t), which measures the 
strength of the coupling. 

Phase locking occurs when |Ao>| < K. The coupling 
has a noticeable effect on the dynamics of Oscillator 1 
only when D\ < O(K), otherwise its dynamics is blurred 
by the internal noise. 

The minimal model for the phase dynamics, Eq. (^), 
can exhibit three different types of dynamical behavior: 

• If there is no phase locking, 0i (t) performs a ran- 
dom walk with a mean drift {dipi/dt) t ~ ujit, see 
Fig. |a. 

• If the phase is locked, it can either be hard-locked, 
i.e. it fluctuates around {d4>\/dt)t ~ u 2 t, see 
Fig. |b, 

• or it is soft-locked, i.e. the fluctuation of 4>\{t) 
around (d(j)i/dt) t ~ io 2 t is intercepted by rapid in- 
creases or decreases of 27r, see Fig. 0c. These jumps 
are called phase slips. 

These different behaviors only occur for D 2 -C D\ which 
is assumed in the following. Otherwise, 0i(i) always per- 
forms a random walk. 

Denoting the average number of phase slips per unit 
time by r, the regime for soft phase locking is, expressed 
in the characteristic time scales, 

D 2 ^r^K , (4) 

whereas for the regime of hard phase locking it is 



Based on these properties of the minimal model, the 
test for phase locking is performed in two steps: A first 
test checks for the frequent occurrence of phase slips. 
If the outcome of this test is negative, a second test is 
performed to distinguish the case of hard phase locking 
from the absence of synchronization. 

In order to perform the first test, the phase increment 

ip T (t)-=Mt + r)-Mt) (6) 

over time intervals of length t is considered. First, an 
estimate of the probability distribution function p{ij) To ) 
for some fixed delay To is determined. Phase slips in 
0i {t) generally lead to a multimodal distribution p(ip To ) 
with maxima separated by 2-7T as displayed in Fig. 0f, if 
To is chosen to be longer than the duration of a phase 
slip 0{K~ V ) J2(| and shorter than the coherence time of 
Oscillator 2: 

K- 1 < r < D^ 1 . (7) 

The estimate of p{4>t ) can be rather noisy. In order 
to obtain a smooth estimator that conserves the possi- 
ble multimodality of p(tp To ), the auto-covariance function 
u([3) oip{ip To ) is considered: 

/oo 
p(VvoMVvo+/?)# , (8) 
-oo 

Multimodality in p{tp To ) translates to multimodality in 
u(0), see Fig. 0g, but the maxima are now located at the 
positions j3 = 2nk (k G Z). 

The degree to which u((3) has such a multimodal struc- 
ture can be quantified by the raw phase-slip index 

^ :v 

Vraw := V u(2/c7t) - U {{2k - 1)tt) . (9) 

v ; fe=-oo 

In the case of soft phase locking, determined by the rela- 
tion between the characteristic time scales of the minimal 
model given in Eq. (0), u((3) can be approximated by a 
sum of Gaussians 

»(«=vfci:«.(x.)exp(-<t|^) , (10) 

where a;(To) gives the probability that in the time in- 
terval To a net number of I phase slips occur, and a is 
the width of the Gaussians. If the internal noise is weak, 
i.e. D\ <C K, these Gaussians are distinct and the contri- 
bution of the minima u((2k — in Eq. (^|) can be ne- 
glected. With J2i a i( T o) — 1 this leads to an upper bound 
Vmax — ao( T o) _1 for Vraw Assuming a Poisson distribu- 
tion for the phase jumps, a[(r ) can be derived analyti- 
cally in dependence of the average number of phase slips 
per unit time r 



r < 0{D 2 ) 



(5) 



ai(r ) = exp(-2rr ) Ii(2rr ) 



(11) 
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FIG. 2: Simulation results, (a) Observed frequency of driven 
oscillator in dependence on frequency of driving oscillator, 
(b) Phase-slip index 77 in dependence on driving frequency 
u)2- (c) Variance of var?/v vs. r. The saturation of vaxtp T 
indicates hard phase locking at u)2 = 1.0 (solid) compared to 
the uncoupled case (dashed). 



FIG. 3: Experimental setup: A pulsator generates a flow pul- 
sating with frequency f pu is- Downstream vortices are formed 
by a shedder bar. The frequency of the vortices f VO rt is de- 
termined via a piezoelectric sensor. 



with /;(.) denoting the modified Bessel function of the 
first kind (22) . The time delay To of the phase increment 
i/j To has to be chosen large enough that a noticeable num- 
ber of phase slips occurs. At the same time it has to 
be small enough that the multimodality of u(/3) is not 
blurred by the phase diffusion of the driving oscillator. 
Numerical simulations show that choosing tq such that 
the Gaussian peak at the origin contributes w 75% to 
u(P) is a good choice. With ao(ro) = 0.75, it follows 
r]max = 1.33. Analytical calculations |^lj reveal that in 
the case of soft locking the value of var(/3), with (3 dis- 
tributed as u{(3), can be estimated as 87r 2 rro. Thus from 
fixing a;(ro) we obtain 

var(/3) = 12.2 . (12) 

Since the actual value of r is not known, Eq. (|lj) is used 
as the condition for choosing t$. 

By definition the lowest conceivable value for rj raw is 
obtained when p(tp T0 ) and consequently u{f3) are single 
Gaussians. With the fixed variance of u((3) this gives a 
minimal value r\ m i n — 0.014 pl[ | for rj raw and a normal- 
ized phase-slip index n can be defined as 



rjrnax ^]min 

which takes values between and 1 depending on how 
well pronounced phase slips are. Large values indicate 
soft phase locking. In the case of hard phase locking, no 
To can be determined such that Eq. (|I^) holds, and n 
cannot be calculated. 

If the test for soft phase locking is negative, a sec- 
ond test has to be performed. Absence of phase slips 
can have two reasons: Either oscillator 1 is not syn- 
chronized to oscillator 2 at all or hard phase locking is 



present. As a means to decide between the two alterna- 
tives, the dependence of var tp T on r is considered. For 
hard phase locking, the phase 4>i(t) is described by a 
stochastic process trapped in a local minimum of the 2tt- 
periodic interaction potential of the two oscillators [B3[ . 
In the simplest case this is an Ornstein-Uhlenbeck pro- 
cess H with cov[0i(t),c/)i(t + r)] = (D/K)exp(-rK), 
where D = D\ + D2 ~ D\. It follows that two regions in 
the dependence of var ip r on r can be distinguished. For 
small t there is an approximately linear increase with 
v&ripT 1=3 2D\T. This saturates at r w K~ x and, for 
larger r, var?/;,- = 2D\/K is constant. Such a saturation 
is the signature of hard phase locking. If oscillator 1 is 
not synchronized at all, 4>\(t) performs a random walk, 
and the variance of tj) T grows linearly in r with slope 2D± 

©■ 

When the test for soft phase locking is applied to the 
general n : m case, the phase <pi(t) obtained from the 
signal has to be multiplied by m prior to computing ip To . 
In the test for hard phase locking, no modifications are 
required. 

The proposed procedure is exemplified by a simula- 
tion study based on the minimal model, Eq. (|^), with 
parameters K = 0.1, D\ — 0.02, D2 = 0, u>i = 1, and 
UJ2 varied between 0.65 and 1.35. Fig. ||a shows the ob- 
served frequency (d4>i/dt) in dependence on u>2, Fig. 
the normalized phase-slip index 77 in dependence on LO2, 
calculated from samples of <fii(t) of length N = 10 6 with 

10 samples per unit time. For UJ2 between 0.96 and 1.04, 

11 cannot be calculated since the condition given in Eq. 
( |l2| ) cannot be fulfilled. Fig. |^c displays the dependence 
of the variance of ip T on r. For the hard phase locked case 
at u>2 = 1) saturates for large values of t, while for 
the uncoupled case with K = 0, the variance increases 
linearly. 
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fusion, while no phase locking is indicated by divergent 
phase diffusion. Both methods were demonstrated on 
experimental data from vortex flow meters. 

The authors thank F. Buhl and P. Riegler for providing 
the flow-meter data and the German BMBF for generous 
support (grant 13N7955). 



FIG. 4: Experimental results, (a) Vortex frequency and (b) 
phase-slip index r\ in dependence on pulsation frequency / pu ia- 
(c) The saturation of vart/v vs. t within the region of hard 
phase locking at a larger pulsation amplitude. 



We apply the proposed procedure to detect phase 
locking in vortex flow meters that are applied to non- 
invasively measure flow velocity of fluids [g5| . These de- 
vices make use of a von Karman vortex street that forms 
behind a shedder bar inserted orthogonal to the flow in 
a pipe. The velocity of the fluid can be determined from 
the frequency of vortex-pair formation / VO rt- In the de- 
vice used here, a piezoelectric sensor inserted downstream 
behind the shedder bar is used to detect the pressure os- 
cillations generated from vortices passing by. A common 
problem of vortex flowmetering is phase locking of the 
vortices to pulsatile flow. 

In the experiments, pulsations with frequency f pu i s 
were generated by adding a periodically blocked flow to a 
steady flow. The degree of flow modulation could be ad- 
justed by the relative contribution of the two, see Fig. ||. 
In a first experiment, the peak-to-peak flow modulation 
was sa 10% of the average flow rate, and the frequency of 
the flow pulsation / pu is was varied. Due to the symmetry 
of the setup |Tq |, the strongest phase locking occurs at 
a frequency ratio of / pu i s : /vort =2:1. In Fig. [|a, the 
shift in / V ort due to phase locking at / pu is ~ 2 / vor t is 
clearly visible. 

The normalized phase-slip index r\, Eq. (|l3|), was com- 
puted from 65 s time series of the sensor signal. As shown 
in Fig. Qb, the presence of phase locking is indicated by a 
pronounced rise in r\. Phase locking was always soft. In 
a second experiment, pulsation was increased to « 40%. 
The condition to calculate 77 were not fulfilled. Figure [|c 
displays the varrp T in dependence on r. The saturation 
for large values of t indicates hard phase locking. 

In summary, based on a minimal model for the phase 
dynamics of unidirectionally coupled oscillators, we pro- 
posed a two-step procedure that provides sensitive and 
specific indicators for soft, hard, and no phase locking. 
Soft phase locking is detected by the presence of phase 
slips, hard phase locking by a suppression of phase dif- 
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